Interplay of rare earth and iron magnetism in i2FeAsO with R = La, Ce, Pr, and Sm: 
A muon spin relaxation study and symmetry analysis 
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We report zero field muon spin relaxation (pSR) measurements on 7?FeAsO with R = La, Ce, 
Pr, and Sm. We study the interaction of the FeAs and R (rare earth) electronic systems in the 
non superconducting magnetically ordered parent compounds of .RFeAsOi-^Fa, superconductors 
via a detailed comparison of the local hyperfine fields at the muon site with available Mossbauer 
spectroscopy and neutron scattering data. These studies provide microscopic evidence of long range 
commensurate magnetic Fe order with the Fe moments not varying by more than 15 % within the 
series 7?FeAsO with R = La, Ce, Pr, and Sm. At low temperatures, long range R magnetic order 
is also observed. Different combined Fe and R magnetic structures are proposed for all compounds 
using the muon site in the crystal structure obtained by electronic potential calculations. Our data 
point to a strong effect of R order on the iron subsystem in the case of different symmetry of Fe 
and R order parameters resulting in a Fe spin reorientation in the R ordered phase in PrFeAsO. 
Our symmetry analysis proves the absence of collinear Fe-7? Heisenberg interactions in i?FeAsO. 
A strong Fe-Ce coupling due to non-Heisenberg anisotropic exchange is found in CeFeAsO which 
results in a large staggered Ce magnetization induced by the magnetically ordered Fe sublattice far 
above T§ e . Finally, we argue that the magnetic 7?-Fe interaction is probably not crucial for the 
observed enhanced superconductivity in .RFeAsOi-zFa, with a magnetic R ion. 

PACS numbers: 75.30.Fv, 74.70.-b, 76.75.+i, 76.80.+y 



I. INTRODUCTION 

The recent discovery of high temperature supercon- 
ductivity in LaFeAsOi-^F^ by Kamihara and coworkers 
has triggered intense research in the Fe-pnictides.— In 
its wake, superconductivity with transition temperatures 
that exceed 50 K have been found in the oxopnictide ma- 
terials in which La is substituted by R = Sm, Ce, Nd, 
Pr, and Gd, respectively^i 3 -^^ Besides the high criti- 
cal temperatures, striking similarities to the properties of 
the high-T c cuprates have been pointed out. As with the 
cuprates, the Fe-pnictides have a layered crystal struc- 
ture with alternating FeAs and RO sheets, where the Fe 
ions are arranged on a simple square lattice X Supercon- 
ductivity emerges in the pnictides when doping the an- 
tifcrromagnctic parent compound either with electrons 
or holes which suppresses the magnetic orderi 14 i 18 These 
similarities raised the hope that cuprates and pnictides 
share a common mechanism for superconductivity, and 
that after 20 year of research on high-T c cuprates the 
Fe-pnictides may provide new insight into the supercon- 
ducting coupling mechanism and verify existing theories 
about high temperature superconductivity. 

In contrast to the cuprates, the non superconducting 
magnetic parent compound is not a Mott-Hubbard in- 
sulator but a metal. Theoretical studies reveal a two- 
dimensional electronic structure with all Fe 3d bands 
contributing to the density of states at the Fermi levels 
Neutron diffraction^ 3 , and local probe techniques like 



Mossbauer spectroscopy and /iS R 12 ' 39 i 40 prove that the 
i?FcAsO parent compounds order in a commensurate 
spin density wave (SDW) magnetic order with a strongly 
reduced ordered Fe moment. Neutron studies suggest 
a columnar magnetic structure with a Fe magnetic mo- 
ment between 0.25 /ib and 0.8 /j,b below ~ 140 K 
that depends on Due to the small size of the ordered 
SDW moment (compared to metallic iron with a moment 
of approx. 2.2 fi-Q per Fe) and the lack of large single 
crystals the temperature dependence of the magnetic or- 
der parameter can be determined by local probe tech- 
niques with a much higher accuracy than with neutron 
scattering^ Note that the magnetic transition in the 
i?FeAsO system is always preceeded by an orthorhombic 
structural distortion which appears at T5, which is about 
10-20 K above T N . 

It is still an open question why the rare earth contain- 
ing systems in the series i^FeAsOi-^F^ with a localized R 
magnetic moment have a higher T c than LaFeAsOi-^F^. 
One suggestion is based on a purely geometrical argu- 
ment. It is argued that the different ionic radii of the 
rare earth elements change the Fe-As-Fe bond angles in 
the Fe-As plane and the FeAs-FeAs intcrplane distance. 
As a consequence, the planar anisotropy of the electronic 
properties may be better for superconductivity in case of 
R = Ce, Pr, Nd, etc.. On the other hand, the electronic 
interaction of the R 4f electrons with the Fe conduction 
band states may be crucial to enhance the density of 
states at the Fermi energy. 
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In this work we examine the interaction of the FcAs 
electronic bands with the rare earth subsystem by a 
detailed comparison of the local hypcrfine fields at the 
muon site and Fe nucleus with neutron scattering results. 
These studies were performed on the undoped magnet- 
ically ordered parent compounds of the i?FeAsOi-. T F x 
superconductors. We report zero field muon spin relax- 
ation measurements on powder samples of i?FcAsO with 
R = La, Ce, Pr, and Sm. We provide microscopic evi- 
dence of static commensurate magnetic order of Fe mo- 
ments. The Neel temperatures and the temperature de- 
pendence of the Fe sublattice magnetization were deter- 
mined with high precision and are compared with avail- 
able Mossbauer spectroscopy and neutron scattering re- 
sults from the literature. In contrast to neutron studies 
our measurements prove that the size of the ordered Fe 
moment is independent of the rare earth ion. For a quan- 
titative analysis of the /iSR spectra the muon site in the 
i?FeAsO crystal structure has been determined by elec- 
tronic potential calculations using a modified Thomas- 
Fermi approach. Spontaneous magnetic ordering of the 
rare earth magnetic moment is observed by /iSR below 
T$ = 4.4(3), 11(1), and 4.66(5) K for R = Cc, Pr, and 
Sm, respectively. Iron and R magnetic structures are pro- 
posed for all compounds on the basis of a detailed sym- 
metry analysis and magnetic dipole field calculations at 
the muon site on the one hand and ^SR, Mossbauer spec- 
troscopy and neutron scattering data on the other. Non- 
collincar magnetic order of Ce and Sm is found in the 
corresponding compounds, which we explain by a weak 
coupling of adjacent R planes in the R-O-R layer. In 
PrFeAsO the /iSR data suggest an Fe spin reorientation 
developing below T^ r , while the Fe order is unaltered 
below T^j in the Sm and Ce compounds. 



In CeFeAsO we find a sizable staggered magnetization 
of the Ce ions induced by the Fe subsystem already far 
above T^ e which amounts to approximately 0.3 hb near 
to Tpj e - We argue that the neglect of the Ce magnetiza- 
tion may have caused the overestimation of the ordered 
Fe moment in recent neutron diffraction studies. Our 
symmetry analysis proves the absence of collinear Fe- 
R Heisenberg interaction in iZFeAsO compounds. Using 
classical and quantum mechanical approaches we deduce 
Fe-Ce and Ce-Ce exchange coupling constants and show 
that the Fe-Ce non-Heisenberg exchange interaction is 
exceptionally strong and of the same order as the Ce-Ce 
exchange interaction. In the Sm and Pr compounds the 
observed coupling between the R and the Fe subsystems 
is found to be much weaker than in CeFeAsO. There- 
fore, we conclude that the magnetic R-Fe interaction is 
probably not crucial for the enhanced superconducting 
transition temperatures in iZFeAsOi-^Fa; with magnetic 
4f ions compared to LaFeAsOi^^Fa;, since only in Ce- 
FeAsO a strong R-Fc coupling is observed. 



II. EXPERIMENTAL 

Polycrystallinc i?FcAsO have been prepared by using 
a two-step solid state reaction method, similar to that 
described by Zhu et al., and annealed in vacuum^ The 
crystal structure and the composition were investigated 
by powder X-ray diffraction. From the X-ray diffraction 
data no impurity phases are inferred. 

In a /iSR experiment nearly 100% spin-polarized 
muons are implanted into the sample one at a time. 
The positively charged /i + thcrmalizc at interstitial lat- 
tice sites, where they act as magnetic micro probes. In 
a magnetic material the muon spin precesses about the 
local magnetic field B at the muon site with the Lar- 
mor frequency / M = 7^/ (2tt)B (muon gyromagnetic ratio 
7 M /27r = 135.5 MHz T" 1 ). With a lifetime of r M = 2.2 ^is 
the muon decays into two neutrinos and a positron, the 
latter being predominantly emitted along the direction of 
the muon spin at the moment of the decay. Measurement 
of both the direction of positron emission as well as the 
time between muon implantation and positron detection 
for an ensemble of several millions of muons provides the 
time evolution of the muon spin polarization P(t) along 
the initial muon spin direction. Magnetic materials with 
commensurate order possess a well-defined local field B 
at the muon site. Therefore, a coherent muon spin pre- 
cession can be observed, which for powder samples has 
the following functional form (see e.g. Ref. @): 
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The occurrence of 2/3 oscillating and 1/3 non-oscillating 
/iSR signal fractions originates from the spatial averag- 
ing in powder samples where 2/3 of the magnetic field 
components are perpendicular to the fi + spin and cause 
a precession, while the 1/3 longitudinal field components 
do not. The relaxation of the oscillation, Ay, is a measure 
of the width of the static field distribution AB = \t/i^ - 
Dynamical effects are also present in At while the relax- 
ation of the second term, A^, is due to dynamic magnetic 
fluctuations only. 

If n magnetically inequivalent muon sites exist in the 
crystallographic or magnetic structure, each of the sites 
contribute to the //SR signal with its weight Pi. In 100% 
magnetically ordered specimens ^2 Pi = 1- Therefore 
/iSR not only provides a highly sensitive measure of the 
magnetic order parameter via internal magnetic fields 
B, but also allows to independently determine the mag- 
netic volume fraction. This is not possible with non-local 
probes such as e.g. neutron diffraction. 



III. MUON SPIN RELAXATION RESULTS 

In Fig. [1] zero field (ZF) /iSR time spectra are shown 
for i?FcAsO with R = La, Ce, Pr, and Sm. At high tem- 
peratures above 150 K no muon spin precession and only 
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a very weak depolarization of the fi + polarization P{t) 
is observed. This weak depolarization and its Gaussian 
shape are typical for a paramagnetic material and reflect 
the occurrence of a small Gaussian-Kubo-Toyabe depo- 
larization originating from the interaction of the fj + spin 
with randomly oriented nuclear magnetic moments.— At 
temperatures below T/v a well-defined spontaneous muon 
spin precession is observed in all compounds indicating a 
well-defined magnetic field at the muon sites. Therefore, 
long range static magnetic order with a commensurate 
magnetic structure is realized in all investigated com- 
pounds of the i?FeAsO series. Accordingly, incommensu- 
rate order or spin glass magnetism can be excluded. Only 
in LaFeAsO, a second /iSR frequency with lower ampli- 
tude P2 ~ 0.30 is observed in addition to the main preces- 
sion signal Pi ?s 0.70. As described above, this indicates 
that two magnetically incquivalent muon stopping sites 
are present in the crystal lattice/magnetic structure. 
The (static) magnetic order develops in the full sample 
volume below T/v as evidenced by the magnetic fj.SK sig- 
nal fraction shown in the inset of Fig. [2] The observed 
5-10% residual non-magnetic signal fraction observed in 
our measurements are due to muons that do not hit the 
sample but stop in the sample holder or cryostat walls. 

The Neel temperatures obtained from the fiSR mea- 
surements (superscript fi) are summarized in table [U 
As we have shown previously for LaFeAsO, the Neel 
temperature T/v of the iron subsystem and the struc- 
tural transition temperature Tg can also be determined 
from anomalies in the temperature dependence of the 
resistivity— For all i?FcAsO compounds investigated 
here, pronounced anomalies have been observed. 10 ' 11 ' 12 ' 15 
The corresponding transition temperatures Tjy and Tg 
are listed also in table |TJ In addition magnetic, of 
the iron and the rare earth (superscript R) subsystems, 

and structural, Tg X , transition temperatures deduced 
from neutron (superscript n) and X-ray scattering (su- 
perscript x) experiments are given for comparison. Note 
that for SmFeAsO no neutron data are available due to 
the high neutron absorption of natural Sm. 

TABLE I: Magnetic and structural transition temperatures of 
iJFeAsO with R = La, Pr, Sm, and Ce. In the cases where no 
reference is given the transition temperatures were measured 
exactly on the samples used in this study . 10 ' 11 ' 12 ' 14 ' 15 All tem- 
peratures are given in Kelvin. See text for the abbreviations 
in the superscripts. 
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Within the experimental error, the Neel temperatures 
determined from fiSK, resistivity and neutron scattering 
agree very well. The structural transition is found to be 
clearly separated by 10-20 K from T/v for all investigated 



compounds. Our highly sensitive /zSR investigations, 
prove the absence of static long range magnetic order 
between T/v and Tg. Also no static magnetic short range 
correlations or disordered magnetism, which would have 
been easily detected by fj.SK, has been observed. How- 
ever, our fjSK data do not rule out a dynamic nematic 
magnetic phase with broken Ising symmetry that has re- 
cently been proposed to develop below Tg ; 20 ' 21 provided 
that the fluctuations are faster than approx. 10 GHz. 



A. Iron magnetic order 

Now we turn our discussion to the temperature depen- 
dence of the Fc sublattice magnetization. We will first 
concentrate on the temperature regime above the static 
order of the rare earth moment. The temperature depen- 
dence of the muon spin precession frequency for iZFeAsO 
with R = La, Ce, Pr, and Sm is shown in Fig. [2] for com- 
parison. All compounds display a very similar tempera- 
ture dependence and absolute value of the fjSK frequency. 
Only the R = Ce compound shows a higher frequency 
and a stronger temperature dependence below T/v. In 
Fig. [3] the yuSR frequency / M (T) is shown together with 
the average magnetic hyperfine field Bhf(T) at the Fe site 
from 57 Fe Mossbauer spectroscopy 2 ^ an d the square root 
of the magnetic Bragg peak intensity yJI(T) of available 
neutron scattering data i 16 ' 17 ' 22 The scale for /^(T) and 
Bhf(T) is the same for all diagrams, and the scales for 
y/l(T) have been adjusted so that the /xSR and neutron 
data coincide as much as possible. In the following, we 
discuss the data for each compound separately. 



1. LaFeAsO 

In LaFeAsO all three observables /^(T), Bhf{T), and 
y/l{T) scale with each other since they all measure the 
size of the ordered Fc moment. A gradual increase of the 
ordered Fe moment observed below T/v indicates a second 
order phase transition. The steep increase of the order 

parameter deviates from the mean-field \J 1 — (T/T/v) 2 
temperature dependence. As we have shown previously, 
this can be understood qualitatively in the framework 
of a four band spin density model ^ The onset of the 
magnetic order is accompanied by a broad static mag- 
netic field distribution AB = At/Tju at the muon site 
(see Sec. [TTJ) that narrows rapidly with decreasing tem- 
perature, as can be seen in Fig. Ufa). The relative width 
of the field distribution A re / = AB /B is proportional to 
(T — T^ e ) _1 , i.e. A re i diverges with T _1 as the tempera- 
ture approaches the Neel temperature of the Fe sublattice 
from below for all systems presented here (not shown). 
Except for SmFeAsO no dynamic magnetic fluctuations 
are detected (Xl = 0) for temperatures below T^ e . 
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FIG. 1: Zero field ^iSR spectra for 7?FeAsO with R — La, Ce, Pr, and Sm for three different temperatures: T > Tat, 
T N >T > Tjy, and T < T§. 
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FIG. 2: Muon spin precession frequency as a function of re- 
duced temperature for 7?FeAsO with 7? = La, Ce, Pr, and 
Sm. Inset: Magnetic signal fraction for 7?FeAsO with R = 
La, Ce, Pr, and Sm. Lines are guides to the eye. 



2. SmFeAsO 

For SmFeAsO no neutron and Mossbauer data are 
available. Qualitatively and quantitatively the tempera- 
ture dependence of the observed /iSR frequency is very 
similar to that of LaFcAsO and consistent with previ- 



ously reported //SR data^i. A sharp increase of the or- 
dered Fe moment is observed below T^ e . Also the satu- 
ration value of /o ~ 23 MHz indicates the same size of the 
ordered Fe moment in LaFcAsO and SmFeAsO assuming 
the same hyperfine coupling constants in these isostruc- 
tural compounds. In contrast to the other systems, mag- 
netic fluctuations in the time window of /iSR are detected 
in this system that cause the dynamic relaxation rate A l 
to increase gradually and to saturate between 10 K and 
30 K. As reported by Khasanov et al., this can be asso- 
ciated with fluctuating Sm magnetic moments due to a 
thermally activated population of Sm crystal electric field 
levels^ The temperature dependence of the dynamic re- 
laxation rate is well described by 



X L (T) X° L CeMEo/k B Ty y ' 

with a saturation value X° L at low temperatures and an 
activation energy Eq that is related to low lying Sm 
crystal electric field levels. The rough agreement of 
the activation energy Eq of the Sm magnetic moment 
fluctuations with the Sm crystal electric field splitting 
has been confirmed by specific heat measurements re- 
ported by Baker and coworkers.— A fit of Eq. (J2]) to 
the longitudinal relaxation rate Al(T) is shown in Fig.Q] 
and yields X° L = 0.172(6) ^s"\ C = 0.001 ms~\ and 
Eq = 44(10) mcV. This value for Eq is approximately two 
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FIG. 3: Main /iSR frequency as a function of reduced temperature T/T^ e together with the average magnetic hyperfine field 
at the Fe site from Mossbauer spectroscopy^ and the square root of the magnetic Bragg peak intensity of available neutron 
scatterin g 16 ' 17 ' 22 data. Note that the scale of f^sR and of Bhf is the same in all graphs. Some axes' labels have been omitted 
for clarity. Typically, error bars of the /j,SK frequency are smaller than the data points. 



times higher than in the oxygen deficient SmFeAsOo.75 
studied by Khasanov et al.i2£ However, one has to take 
into account that the oxygen deficiency in SmFeAsOo.75 
causes a higher defect density and changes of the lattice 
parameters^ Therefore, a different activation energy Eq 
for the undoped SmFeAsO compared to SmFeAsOo.75 is 
possible. 



systematically reduced by a small amount compared to 
Bh / (T) . This phenomenon can be explained by the same 
mechanism found in the Ce system (see below and Sec- 
tion lVip . In the case of PrFeAsO the muon spin polariza- 
tion function (JTJ) did not approximate the data very well 
and a non-zero phase <p of the oscillation and a gener- 
alized exponential relaxation function exp(— (At*)") had 
to be used to describe the data. 



3. PrFeAsO 

Apart from a slightly reduced Neel temperature, 
PrFeAsO shows the same temperature dependence and 
saturation value of the Fe sublattice magnetization, i.e. 
fn(T) as SmFeAsO. Note that again the saturation value 
for /„ and f are very similar to that of LaFeAsO indi- 
cating the same size of the ordered Fe moment. A close 
comparison of f^(T) and Bh/(T) reveals that / M (T) is 



4. CeFeAsO 

A qualitatively different behavior is observed for the 
CeFeAsO compound. Neutron diffraction as well as /xSR 
data do not scale with the observed hyperfine field Bhf 
at the Fe site. The magnetic Bragg intensity as well as 
the internal magnetic field at the muon site measured by 
the muon precession frequency continuously increase be- 
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FIG. 4: top: Transverse relaxation rate A-r(T) for 7?FeAsO 
with R — La, Ce, Pr, and Sm. Lines are guides to the eye. 
bottom: Longitudinal relaxation rate Xl(T) for iJFeAsO 
with 7? = Sm only. The solid line is a fit of equation [2] to 
the data. The dashed line is a guide to the eye 



low Tjv. Only the Mossbauer hypcrfinc field displays the 
same rapid saturation below T/v with the same ordered 
moment as observed for the La and Pr compounds.— 
Therefore we conclude that the fiSK as well as the neu- 
tron data do not solely measure the Fe sublattice magne- 
tization, but also contain a significant contribution from 
the Ce sublattice. The 57 Fe Mossbauer spectroscopy 
provides the most accurate measurement of the on-site 
Fe sublattice magnetization without sizable contribution 
from the rare earth moments due to a weak transferred 
hyperfine coupling. In contrast, a magnetization on the 
rare earth site induced by the Fe subsystem has the same 
symmetry as the Fe order and therefore contributes to the 
same Bragg peaks. In principle, neutron scattering can 
distinguish between the different contributions from the 
Fe and the Ce sublattice by fitting the different magnetic 
form-factors. This has not been done for T > in the 
neutron studie a 16 i 17 ' 22 shown in Fig. O where the whole 
magnetic intensity has been attributed to originate from 
Fe moments alone. Similar to the neutron data, the local 
field at the muon site also contains a contribution from 
the rare earth magnetic sublattice. In the following, we 
model the above mentioned contribution of the induced 
Ce moment to the field at the muon site by an additional 



Curie-Weiss contribution. This can be interpreted as 
a local magnetization of paramagnetic moments of the 
Ce sublattice induced by the molecular field generated 
by the Fe sublattice. In turn, the induced Ce sublattice 
magnetization creates a dipole field at the muon site. The 
temperature dependence of the Fe molecular field has to 
be proportional to the Fe sublattice magnetization, i.e. 
the muon spin precession frequency /^(T) observed in 
the other i?FcAsO systems. This is plausible because 
we have shown that J M (T) above the rare earth ordering 
temperature is almost independent of the rare earth 
ion. Thus we chose the following function to describe our 
data in a first approximation: 



1 - 



T 

rpFe 
1 N 



C 



T 



(3) 



The first term in the last square brackets is used to de- 
scribe the contribution of the Fe sublattice, the second 
term is the additional Curie- Weiss contribution of the Ce 
sublattice to the magnetic field B(T) = 2ir f u(T) / '7 „ at 
the muon site. Here, C describes the hyperfine coupling 
constant of the Ce moments with the muon spin. A fit 
of this function to the fiSR frequency / M (T) obtained for 
CeFeAsO is shown in Fig. [3] and the two contributions 
are highlighted as described in the legend. This simple 
model describes the data reasonably well for tempera- 
tures between 10 K and up to T^ e — 137 K, yielding 
fo = 25.7(5) MHz, a = 2.4(4), (3 = 0.24(1). The con- 
stants C = 2.5(5) K, and 6 = —3(0.5) K were obtained 
by restricting the fit range to the most relevant tempera- 
ture region between 10 K and 50 K. The exponents a and 
(3 are close to the results obtained for the other systems 
without the Curie-Weiss contribution. The enhanced fo 
could be due to a slightly different muon site compared 
to other i?FeAsO compounds. The deviation from this 
behavior below 10 K is attributed to the growth of anti- 
ferromagnetic correlations in the vicinity of T^ e . 

The exceptionally strong coupling of the Ce is already 
reasonable considering the ground state properties of the 
R ion in the crystal electric field (CEF). The ground state 
of the Kramers ion Ce 3+ is the J-multiplet 2 F 5 / 2 - The 
Sm 3+ Kramers ion has 6 H 5 / 2 as lowest J-multiplet with 
J = 5/2. The non-Kramers ion Pr 3+ stays in the 3 i?4 
ground multiplet with J = 4. The magnetic behavior 
of the R ions can be understood qualitatively from the 
susceptibility of free R ions: 



Xo(T) 



g>|j(J+l) 
3T 
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Thus, because of the different g-factors (gj(Ce) = 6/7, 
g.;(Sm) = 2/7, and gj(Pr) = 4/5) one expects at least 
one order of magnitude less induced Sm magnetic mo- 
ment compared to Ce (for equal magnitudes of Fe-Sm 
and Fe-Ce coupling constants, see section |VB|) . Prob- 
ably, due to this feature, we observe very similar ^tSR 
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FIG. 5: Structure and electrostatic potential energy map of the (llO)-plane for a muon in the Cmma orthorhombic phase of 
CeFeAsO. The A and B muon positions are marked by corresponding letters. The potential energy is given in atomic units. 
The origin of the unit cell has been moved for clarity. 



spectra in the Sm, Pr, and La compounds in the tem- 
perature range T$ < T < T^. since the exact sequence 
of CEF levels in PrFeAsO is not known up to now. it is 
more difficult to predict its behavior. 

The above qualitative discussion neglects higher CEF 
levels and the Ce-Ce interactions. A more detailed and 
quantitative discussion of the coupling of the Ce sub- 
lattice to the Fe sublattice using classical and quantum 
approaches will be given in Sec. Mand IVII 



B. Rare earth magnetic order 

For R = Pr and Ce the onset of the rare earth mag- 
netic order 11(1) K and 4.4(3) K, respectively causes a 
second order, i.e. smooth decrease of the muon preces- 
sion frequency ffi(T) by approx. «2.2 MHz at «2 K 
in both cases, as can be seen in Fig. [3J However, in 
PrFeAsO the Pr magnetic order is accompanied by a 
maximum in the width AB of the magnetic field dis- 
tribution, i.e. the relaxation rate Xt (see inset of Fig. [4j 
that rapidly decreases at lower temperatures, whereas in 
CeFeAsO the Ce magnetic order causes the field distribu- 
tion to broaden monotonically. Note that the magnetic 
field at the muon site caused by the static order of the 
Ce sublattice leads to a decrease of the observed fiSR 
frequency while the Fe order induced magnetization of 
the Ce sublattice above T^ e causes an increase of the 
frequency f^{T). Therefore, it is evident that the Ce 
moments order in a different structure than induced by 
the Fe magnetic system. 

SmFeAsO shows first order Sm magnetic order, i.e. the 
Sm magnetic order parameter is discontinuous at T^ m = 
4.66(5) K. Contrary to the Ce and Pr magnetic order 



the Sm magnetic order causes the appearance of two 
satellite muon precession frequencies /^ ,2 (r) = f®(T) ± 
8.0(5) MHz in addition to the main frequency f®(T) ob- 
served above T^ m . The two satellites and the main fre- 
quency have signal fractions of 1.0(3):1.0(3):4.3(6). This 
shows that the Sm magnetic order has different symme- 
try compared to the Fe magnetic order, i.e. it causes a 
change of the magnetic unit cell. A detailed discussion 
of R magnetic structures and its interplay with the Fe 
sublattice will be given in Sec. IVII IVHI and I Villi 



IV. DETERMINATION OF THE MUON SITE 

To determine the contributions from Fe and R mag- 
netic order to the local magnetic field at the muon site 
it is necessary to determine the muon site in the lat- 
tice. For this purpose we used a modified Thomas Fermi 
approach^ and available structural data. The goal is 
to determine a self consistent distribution of the valence 
electron density from which the electrostatic potential 
can be deduced. The local minima of this potential at in- 
terstitial positions are regarded as possible stopping sites 
for muons. We verified the applicability of our approach 
by comparison of numerical results with experimentally 
determined muon sites in iJFeO^^I and by a successful 
interpretation of /xSR spectra of the complex magnetic 
structures in layered cobaltites i?BaCo20s.5i^ 

A potential map of CeFeAsO in the Cmma orthorhom- 
bic phase was calculated using structural data*£ at 1.4 K 
and is shown in Fig. [5] The calculations have been done 
without taking into account the host lattice relaxation 
around the muon. We observe two types of possible muon 
positions which are labeled A and B. The A type posi- 
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tion is located on the line connecting the Ce and As ions 
along the c-direction. Note that a similar location of the 
point with deepest potential was calculated for LaOcAs 
in the tetragonal phase using the general potential lin- 
earized augmented plane-wave method and local density 
approximation.— The A type position with coordinates 
(0,1/4^) has 4g local point symmetry mm2, i.e. the 
same as the R sites. The B type position has a gen- 
eral 16o local point symmetry with (x,y,zs) coordinates. 
For CeFeAsO the coordinates of A are (0,1/4,0.41) and 
(0.44,0.19,0.04) for B. Note that the points B are located 
at oxygen ions which are typical points for muons in many 
oxides. The za and zb coordinates vary in the third dec- 
imal if Ce is replaced by La, Pr or Sm. 



A second important result is the comparison of the 
calculated muon precession frequency, i.e. the local 
magnetic dipole field at the muon site with the exper- 
imentally determined muon precession frequency /o = 
/ M (T — > 0) for LaFeAsO. With the experimentally de- 
termined Fe magnetic momenl!^ 3 - of 0.36 /xb we calculate 
f ,4 = 12.4 MHz for site A and f 0tB = 1.3 MHz for site 
B. As shown by Klauss et al. in Ref. [HI and Section [TTT1 
of this work, in LaFeAsO two muon frequencies arc ob- 
served. One frequency with 23 MHz originating from the 
mayor volume fraction (P = 0.7) and one lower frequency 
with 3 MHz which develops from a strongly damped sig- 
nal below approximately 70 K. In i?FcAsO with R = 
Sm, Pr, and Ce only the high frequency is present in the 
/iSR spectra. We conclude that site A is the main muon 
site since this gives the correct order of magnitude for 
fo- Probably the site B is also partially populated at 
low temperatures in the LaFeAsO. The fact that we ob- 
tain a 46 % smaller value than the experimental result is 
reasonable since our calculation only takes into account 
local dipole fields and neglects contact hypcrfine contri- 
butions. Similar differences are found in La2Cu04, the 
antifcrromagnctic parent compound of the 214-cupratc 
superconductors^ As will be shown below, the magnetic 
field caused by the Fe magnetic order is directed along 
the crystallographic c-axis. This is in agreement with 
recent [iSH experiments on Fe pnictide single crystals^ 
In the following we use only muon site A and a renor- 
malization factor of 1.86 for the local field caused by Fe 
magnetic ordering at this muon site to account for the 
contact hypcrfine field for all R. This assumption is jus- 
tified from our calculation of the electronic charge den- 
sity distribution that shows very similar results for all R. 
This manifests itself in the almost identical position of 
muon site A in compounds with different R. In its turn 
this means that the renormalization due to the contact 
hypcrfine interaction, which depends only on the elec- 
tron density at the muon site, should be nearly the same 
within this series. 



V. SYMMETRY ANALYSIS OF FE AND R 
MAGNETIC ORDER PARAMETERS IN i?FEASO 

A. Translational symmetry and magnetic modes in 
.RFeAsO compounds 

Iron magnetic order in i?FeAsO sets in about 10-20 K 
below the tetragonal to orthorhombic structural phase 
transition^ The space group of the paramagnetic phase 
is Cmma with Fe and R ions occupying 4b and and 4g 
position respectively. Neutron studies of i?FeAsO com- 
pounds revealed numerous magnetic Bragg peaks all of 
which can be reduced just to three types of magnetic 
propagation vectors i 11 ^ 3 ' 16 ^ 7 ' 22 ^ 3 ' 47 These vectors in an 
orthorhombic setting are: 

K 7 = (1,0,1/2) 

K n = (1,0,0) (5) 
K IH = (0,0,1/2). 

Here we use an orthorhombic primitive cell with the unit 
cell vectors ao = 2r x e a , bo = 2r y eb, and Co = 2r z e c 
where t x , t v , and 2t z are the distances between nearest 
neighbor Fe ions along a, b and c directions, denoted by 
the unit vectors e, , respectively. Note that in this setting 
nuclear Bragg peaks of the Cmma spacegroup are either 
(2H,2K,L) or (2H + 1,2K + 1,L). Also, a 'magnetic 
setting' with a m = ao, b m = bo, and c m = 2co is often 
used. 

The primitive cells of the magnetic structures induced 
by each of the K 7 are not orthorhombic. Each of the 
magnetic primitive cells contains four magnetically in- 
equivalent 7?-, Fe-, and muon-sites. To analyze the sym- 
metry of the magnetic structures in a unified way we will 
label these four positions with numbers a = 1,2,3,4, 
with corresponding ordered magnetic moments m a for 
the R- and Fe-positions. To identify the different posi- 
tions a in the magnetic cell with positions (3 = 1, 2, . . . , 8 
in the nuclear cell (doubled along the c-direction) we use 
the mapping {/3i,/32} > ol where a is one of the four 
possible positions in the magnetic primitive cell and fti a 
position in the nuclear cell as shown in Fig. [SJ 

The Kj propagation vector induces a magnetic lattice 
of I c mmm type in which the primitive cell is built on 
the translations ai = ao, &2 = bo, a,3 = 1/2 (ao + 
bo) + c o an d the anti-translation along the z-axis a' = 
co^ Thus, in Kj-typc magnetic structures the following 
positions are equivalent and mapped as denoted by W: 
{1, 7} i-> 1, {2, 8} i-> 2, {3, 5} i-> 3, and {4, 6} i-> 4. 

The primitive cell of the Pcmmm magnetic lattice in- 
duced by the K// propagation vector is constructed from 
the translations ai = ao, &2 = bo, a3 = Co and the 
anti translation a' = 1/2 (ao + bo)^ Then we obtain 
the identity of the positions: {1,5} i— > 1, {2,6} i— ► 2, 
{3,7} h-> 3, and {4,8} i-> 4. 

Finally, for the 'Km propagation vector the magnetic 
lattice is C c mmm and the magnetic primitive cell can 
be built from the translations: ai = l/2(ao + bo), 
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FIG. 6: Enumeration of positions used to describe different 
types of Fe and R magnetic order, (a) R and Fe positions and 
b) muon site A. 



&2 = 1/2 (ao — bo), a3 = 2co>^ Therefore we have 
the mappings: {1,3} i— > 1, {2,4} i— » 2, {5,7} i— > 3, and 
{6,8} i-> 4. 

In order to investigate the magnetic interactions in 
different R and Fe magnetic structures with the trans- 
lational symmetry of the propagation vector K 7 (7 € 
{0,1,11,111}) it is necessary to calculate the magnetic 
modes, i.e. basis functions of the irreducible represen- 
tations (IR) of the propagation vector little groups G 7 . 
These groups are the same for all propagation vectors 
([5|) and they are also identical with the groups for ferro- 
magnetic order, i.e. propagation vector Ko = (0,0,0).— 
Following a method described by Bertout, and Izyumov 
and Naish we introduce possible magnetic modes F and 
hi, i=l,2,3 as linear combinations of Fe (or R) sublat- 
tice magnetic moments m Q , where a denotes a particu- 
lar atom (see above) in the respective magnetic primitive 
ce u« 

F = l/4(mi + m2 + m3 + 1114) 
Li = l/4(mi + m2 — m3 — 1114) 
L2 = l/4(mi — m2 + m3 — 1114) 
L3 = l/4(mi — m2 — m3 + m4). (6) 

The components of these linear combinations are the ba- 
sis functions of the eight one-dimensional IR t„ with 
1 < v < 8 for the G group (see table 32 in Ref. [13). 
The results of the symmetry analysis for the magnetic 
moments located at the 4b positions (Fe ions) and 4g 
positions (i? ions and muon sites) are shown in Table HT1 



In accordance with neutron data the iron magnetic 
subsystem orders antifcrromagnctically along the in- 
direction and fcrromagnctically along the ^-direction 
with doubling of the unit cell along the c-axis for the R 
= La compound (magnetic propagation vector K/) and 
without doubling for R = Ce and Pr (magnetic propa- 
gation vector K77).— In all these structures Fe magnetic 
moments are directed along the a-axis. These magnetic 
structures are described by L 3x (Kj) and L 3x (Ku) non- 
zero magnetic order parameters shown in Table [III 

The onset of a given type of magnetic order lowers 
the symmetry of the paramagnetic phase and creates a 
magnetic symmetry. The orientation of the magnetic mo- 
ments and the distribution of local magnetic fields in the 
magnetic cell have the same symmetry. To determine the 
orientation of the magnetic field at a particular muon site 
we assign an artificial magnetic degree of freedom to this 
site. The set of magnetic degrees of freedom for different 
points (Wyckoff positions) forms a magnetic representa- 
tion. A standard decomposition of this magnetic rep- 
resentation into irreducible representations of the space 
group allows the analysis of the symmetry of magnetic 
field distributions at the muon sites. This symmetry 
must be compatible with the given space group and the 
same as the symmetry of the magnetic order parameter. 
In other words, the symmetry of the magnetic field dis- 
tribution at the muon site must belong to the same IR 
as the magnetic order parameter. 

From Table [III it becomes self-evident that the 
L 3x (K.u)~type of magnetic structure of iron moments 
creates L§ z (Kjj)-type molecular fields at the rare earth 
and A-type muon sites. Therefore, a magnetic field dis- 
tribution of any symmetry that induces a staggcred-type 
ordering of the R magnetic moments leads to a mag- 
netic field distribution at the A-type muon site of same 
symmetry. Due to symmetry reasons, some types of Fe 
magnetic order, for instance L 3y {K.i), do not create mag- 
netic fields at the 4g positions. Thus, the presence of 
/iSR signals in the Fe antiferromagnetically (AFM) or- 
dered phase of R = La, Pr, Ce and Sm compounds ex- 
cludes the possibility for AFM magnetic order along the 
a-axis and ferromagnetic (FM) order along &-axis with 
iron magnetic moments directed along 6-axis. 

The L 3x (Ki)~type of Fe order and the corresponding- 
direction of the magnetic field at the A-type muon site 
are shown in Fig. [7] The dipole fields at the A-type 
muon position for all types of Fe and R magnetic order 
are given in the Appendix A. 



B. Symmetry of the Fe— R magnetic interactions 

To start with the strongest exchange interaction we 
analyze the symmetry of permutation modes ([3]) or ex- 
change multiplets for Fe and R magnetic subsystems.— 
The permutation modes for the different translational 
symmetries K 7 are listed together with the correspond- 
ing irreducible representations in table fTTTl 
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TABLE II: Magnetic modes of the magnetic propagation vectors K 7 , i.e. basis functions of the irreducible representations t v 
(see text). The indexes indicate the non-zero Cartesian components of the vectors ([6]), i.e. for example L\ z is the z-component 
of Lr. 
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FIG. 7: (a) L3a;(Kj)-type Fe magnetic order realized in 
LaFeAsO and probably in SmFeAsO and its magnetic dipole 
field direction at the muon site A. (b) Ls^Kr/) order observed 
in CeFeAsO and PrFeAsO. 



TABLE III: Permutation modes, i.e. exchange multiplets of 
the four magnetic propagation vectors. 
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Ko and Hm translational symmetry. The respective 
Fc-R magnetic exchange Hamiltonian has the form: 

Ulr R = ... + J Fe ^(K„)(F.F fl ) 

+J 1 Fe - fi (K /// )(L 1 (K /// )Lf(K /// )) (7) 

However, for the cases of K/ and K// translational sym- 
metry the Fe and R subsystems can interact by the non- 
Heisenberg exchange (see Table ITTT|) . The part of these 
anisotropic Fe-R interactions relevant for the following 
considerations is given below: 

H F a *Z* = ... + /3 F ;- fl (K J )L3.(K / )L 3 «(K / ) + 
+I^~ R (K I )L 3z (K I )L R (K I ) + 
+4 x r R (Kn)L 3x (K II )L£(K II ) + 
+I^- R (K II )L 3z (K II )L R x (K II ) + 
+ ... (8) 

Consequently, the onset of the -^3 K (Kj/jj) type of Fe 
magnetic order creates an effective staggered magnetic 
field at the R site along the z-direction. The magni- 
tude of this field is proportional to the value of the Fe-R 
coupling constant I 3xz ~ R and the iron subsystem order 
parameter L 3x determines its temperature dependence. 
In the following we will show that this field exceeds a re- 
spective dipole field at the R sites by at least one order of 
magnitude. The mutual orientation of the L 3x and L R Z 
vectors depends on the sign of the Fe-R coupling con- 
stant I 3x e ~ . The orientation of the induced magnetic 
moment on the R site is shown in Fig. [5] 

The magnitude of the z-component of the rare earth 
magnetic moments m R z (a = 1 — 4) is determined by the 
exchange field and the R ion crystal electric field (CEF). 



Since there arc no permutation modes of the Fe- and 
-R-subsystems which belong to the same irreducible rep- 
resentation there are no Heisenberg exchange interactions 
between the R- and Fe-subsystems for the case of K/ or 
K.jj translational symmetry of the Fe magnetic order. 
This interaction exists only for magnetic structures with 



VI. DETERMINATION OF THE MAGNETIC 
STRUCTURES AND THE FE-R COUPLING 
CONSTANTS FROM THE TEMPERATURE 
DEPENDENCE OF THE /LtSR RESPONSE 

Below the Neel temperature the /xSR response of 
i?FcAsO is mainly determined by magnetic dipole fields 
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FIG. 8: Direction of the magnetic moment on the R site 
induced by the \jz x {K-ii) Fe order in a) CeFeAsO and b) 
PrFeAsO. The sign of the coupling constant I^~ R is dif- 
ferent in the Ce and the Pr compound resulting in opposite 
orientation of the induced moments. 



created by the Fe and R subsystems on the A-type muon 
site. It can be calculated by: 



B Z {A) = 2.58 •L 3a (K / ) + 1.45 
B Z (A) = 2.57.L to (K/j) + l. 



R 

3z 
E 
?,z 



(Ki 



(9) 



If the Ljj(K 7 ) arc given in units of /ib (sec equation ((6])) 
the magnetic dipole field B Z (A) is obtained in units of 
0.1 T. These equations demonstrate a surprisingly close 
agreement of dipole fields for Fe magnetic order with 
(K/) and without (Kn) doubling along z-axis. Note 
that the constants in equations @ only vary in the third 
decimal if Ce is replaced by La, Pr or Sm. This explains 
the similar /iSR spectra of LaFeAsO with L^ x {K.i)— , and 
CeFeAsO and PrFeAsO both with L 3a; (Kjj)-type of Fe 
antiferromagnetic order. The high temperature fit of the 
/iSR frequency gives roughly the same saturation value 
fo = 7;,/2ttB = 23.0(5) MHz. This corresponds to a 
value of 0.36 /xb for the Fe saturation moment if one 
takes into account contact hypcrfme fields by applying 
the renormalization factor of 1.86 as in LaFeAsO. Be- 
cause of the similar order parameters of the Fe mag- 
netic order found in SmFeAsO, LaFeAsO, and PrFeAsO 
we conclude that a Z^-type of Fe antiferromagnetic or- 
der is also realized in SmFeAsO. However, /xSR studies 
alone cannot distinguish between the possible transla- 
tional symmetries, i.e. K/ and K// of the order param- 
eter. 



1. CeFeAsO 

The /xSR response drastically changes below the Ce 
ordering temperature T^ e . Contrary to the naive expec- 
tation that the Cc magnetic order would increase the lo- 
cal magnetic field at the muon site and therefore the /iSR 
frequency our experimental results show a decrease of the 
/iSR frequency below T^ e . According to our analysis of 
the magnetic field distribution at the A-type muon site 
this behavior indicates a breaking of the L^ z symmetry of 
the induced Ce order. Neutron diffraction studies of Ce- 
FeAsO revealed the appearance of K///-type magnetic 
Bragg peaks below T^ e in addition to K//-type peaks 
also present above Tpj e J£ The Ce magnetic structure pro- 
posed by Zhao et al. is a non-collinear arrangement of 
Ce moments in the afe-plane^i This structure is shown in 
FigJS] a) and can be described as a linear combination of 
two order parameters: L^ x (Kui) + L§ y (Ku). This new 
type of R order produces the following magnetic field at 
the A-type muon site for temperatures below T^ e : 



B X (A) 

By (A) 

B Z (A) 



-0.40 • L%(Km) 



= -0.88-.Lgf(K„) 



2.57- 



fo 



13.55 



lM-Lg B (K ZI ) (10) 



Here /o=25.7 MHz accounts for the contribution of the 
iron subsystem, which, we suppose, preserves its mag- 
netic structure found above T^ e . Despite different trans- 
lational symmetry of the Ce order parameters, the mag- 
netic fields at all eight muon sites have the same mod- 
ulus p0[) . therefore no additional frequencies appear at 
temperatures below T^ e . With the Ce magnetic mo- 
ment of 0.83 /Ub at 1.7 K reported by Zhao et alJ£ and 
supposing L£f=0 we then obtain a /xSR frequency of 
/ = 26.9 MHz instead of the experimentally observed 
29.8 MHz at 1.9 K. Also the experiment reveals a remark- 
able decrease of the observed /xSR frequency compared 
to maximal value of 32 MHz at 4 K. The decrease of 
the /xSR frequency is not as rapid as expected, it devel- 
ops only gradually below T^ e . As follows from equation 
(JTDJ), the magnetic dipole fields at the muon sites cre- 
ated by the Ce magnetic moments lie in the a6-plane 
and are rather weak compared to the dipole field cre- 
ated by the Fe moments. The slow decrease of the ^SR 
frequency below Tjy therefore reflects the gradual dis- 
appearance of the Fe induced L£f type magnetization of 
the Ce moments in the AFM Ce phase below T^ e . One 
can conclude that the Ce-Ce coupling constant which 
creates the L^(K nl ) + L^(K n ) AFM Ce order pa- 
rameter cannot immediately suppress the induced L^ z 
magnetization. As will be shown below the Ce-Ce cou- 
pling constant J^(Ku) is of the same order as the Fe- 
Ce coupling constant I^~ Ce (Kjj). Using equation (fl7J| 
we estimate the remanent order parameter L£, e =0.14/ZB 
at 1.9 K. 

Due to the significant contribution of Cc dipole fields, 
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(a) CeFeAsO (Ref. 16) 




(b) PrFeAsO 




(c) SmFeAsO 



'-+4. 



j 




FIG. 9: (a) Collinear Fe magnetic order (Kn) and non-collinear Ce order in CeFeAsO as proposed by Zhao et al.— which is in 
accordance with the fiSH measurements, (b) Magnetic structure of the Pr and Fe subsystems with K// translation symmetry 
which satisfactorily describes /iSR and Mossbauer data. Non-zero are the L^ x cos 6 + L- iz sin 9 of Fe and L% x of Pr order 
parameters. The Fe and Pr subsystems remain in the L3 and type of exchange magnetic order, respectively, (c) The 
proposed Sm structure which explains the observed /iSR spectra. Fe magnetic ordering has K/ translational symmetry with 
a doubling along the c-axis, while the Sm magnetic order possesses K.n translational symmetry without a doubling along the 
c-axis) . 



the [iSR frequency fp,(T) is unusually high already far 
above T^ e in CeFeAsO. This has been discussed phe- 
nomenologically in Sec. IIIII Here we carry out a more 
elaborate analysis of this observation. The lowest crys- 
tal electric field (CEF) excitations of CeFeAsO and Sm- 
FeAsO consist of three Kramers doublets. With the onset 
of antiferromagnetic iron order these doublets split. In- 
elastic neutron scattering studies of CeFeAsO revealed 
that in the paramagnetic phase the energies of the first 
and the second excited doublets are 216 K and 785 K 
respectively^ The splitting at T = 7 K amounts to 
10.8 K for the ground state doublet and 37.1 K and 66.1 K 
for the first excited levels. Using these experimentally 
deduced parameters and the wave functions of the cor- 
responding CEF levels allows us to determine the Fe-Ce 
and Ce-Ce exchange interaction coupling constants. In 
the Fe antiferromagnetic phase above T^ e , they are ob- 
tained by fitting the theoretical temperature dependence 
of the Ce magnetic moment m^ e (T) to the experimental 
data. In Fig. [10] the experimental value of the Ce single 
ion magnetic moment m^ e (T) is shown. It is calculated 
from the fiSR frequency f exp with the help of equation 
© after subtracting the Fe contribution /o (T) (see inset 
of Fig.[10]): 

Lg e (T) = m c z e {T) = [/ eap (T)-/ (T)]/13.55-1.48 (11) 

In this effective field approximation the Fe-Ce coupling 
constant I 3x e ~ (Ku) and the Ce-Ce coupling constant 
J^ z e (Kjj) both create a staggered effective field B z '* {T) 



on the Ce sites. Its interaction with the Ce magnetic 
moments lifts the degeneracy of the Kramers doublets. 
The corresponding effective field on the Ce site is given 
by: 

B e Jf(T) = B^~ Ce (T) + J 3z m^(T) (12) 

where Bf e " Ce (T) = li x e - Ce {K n )L Zx {T) and J 3z = 
2J§ Z (K U ) with / and J in units of T/^b- We neglect 
contributions from Fe-Ce and Ce-Ce dipolc fields to the 
splitting of Kramers doublets under the assumption that 
the Fe-Ce and Ce-Ce (non-)Heisenberg exchange fields 
(see equation IT2"j) are much stronger. The Fe-Ce stag- 
gered field Bf e ~ Ce has to have the same temperature 
dependence as fa{T), i.e. the averaged value of the iron 
order parameter L 3x (T). Therefore, it can be modeled 
as Bf e - Ce (T) = B ■ [1 - {T/T c ) a f with a = 2.4 and 
(3 = 0.24 as determined in Sec. IIIII 

According to the analysis of inelastic neutron scatter- 
ing data££ the lowest Kramers doublet possesses the or- 
bital momentum J z = ±1/2 and does not contain any 
admixtures of the J z = ±3/2, ±5/2 states. One can esti- 
mate the magnitude of the effective field from the split- 
ting of the ground Kramers doublet at T = 7 K by using 
the relation: 

g Jf i B Bff(T) = k B A (13) 

where g,j = 6/7 is the g- factor of the free Ce 3+ ion and 
A = 10.8 K the splitting of the ground state doublet. If 
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we account for just the Fe-Ce interaction it yields Bq = 
18.75 T. This large Fe-Cc effective field indicates that 
dipole-dipole interactions (see Appendix B) do not play a 
significant role for the magnetization of the Ce sublattice 
by the ordered Fe sublattice above Tjy e . 

The temperature dependence of the Ce moment in the 
effective field has been calculated in a quantum mechan- 
ical approach according to equation (10) in Ref. |42| by 
calculating the thermal population of all experimentally 
determined CEF levels^ The result is shown in Fig. |TD] 
as a green line. From our calculation it is clear that there 
is a deviation from Curie- Weiss behavior due to a sizable 
contribution to m^ e from higher Kramers doublets for 
temperatures above 50 K. In Sec.|TTT]we did a Curie- Weiss 
approximation of m^ e by using Eq. ([3]). A Curie- Weiss 
description neglects higher crystal field levels. This is the 
reason why we restricted the fit to the temperature range 
between 10 and 50 K in the first approximation done in 

Sec. In] 

On the other hand, a systematic deviation between 
m^ e obtained from the quantum approach and the ex- 
perimental data is observed in the low temperature re- 
gion between 10 and 50 K. This indicates that the Ce-Ce 
interaction, which is neglected by this approach is signif- 
icant in this temperature region and has to be taken into 
account for a proper description of the data. Therefore, 
we fit the Ce magnetic moment with contributions from 
the Fe-Ce and Ce-Ce exchange interactions by using the 
usual Brillouin function. In the case of S=l/2 this leads 
to the following equation: 

mf e (T) = ^tanh(|g^/(T)) (14) 

where B & J^ is given in Eq. (TT!?|) which should fulfill the 
boundary condition given by Eq. (|13|) . This approach 
neglects the higher crystal field levels and therefore can 
only be applied at temperatures below 50 K. From a fit 
in the low temperature region we obtain the microscopic 
parameters B = 26.8(5) T and J 3z = -24.3(5) T//z s . 
The fit is shown as a magenta line in the Fig. [TDJ Us- 
ing 0.41 \ib for iron saturated moment one can estimate 
a rather large Fe-Ce coupling constant I%£ z (Ku) = 
65.3 T/fiB- The negative sign of the Ce-Ce coupling 
constant J^ z e (Kjj) = —12.15 T/hb indicates that the 
type of Ce exchange magnetic order is energy fa- 
vor able. From T$=4.4 K one can estimate the value of 
the J^ y e (Kji) coupling constant. It is responsible for the 
appearance of the in-plane L^ y e order parameter^ This 
rough estimate gives J^y(Ku) = —14.9 T/^ib- Note, 
that the big difference between Ce-Ce exchange coupling 
constants along z- and y-directions points to a strong 
easy-plane anisotropy in the Ce subsystem. 

In conclusion, for CeFeAsO we found a sizable stag- 
gered magnetization of the Ce ions induced by the Fe 
subsystem already far above T^ e which amounts to ap- 
proximately 0.3 near to Tjy e . With the help of a 
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FIG. 10: Temperature dependence of the induced Ce mag- 
netic moment mf e (T) for T > T§ e . m^ e (T) has been ex- 
tracted from the fiSH frequency by subtracting the contribu- 
tion from the Fe subsystem (inset). The data are fitted by a 
quantum and a Brillouin approach valid in different temper- 
ature regions, see text. 



symmetry analysis as well as classical and quantum me- 
chanical approaches we were able to deduce Fe-Ce and 
Ce-Ce exchange coupling constants. We demonstrated 
that in the Ce compound the Fe-Ce non-Hcisenberg ex- 
change interaction is exceptionally strong and of the same 
order as the Ce-Ce exchange interaction. The observed 
//SR frequency in the Ce ordered phase is in agreement 
with the magnetic structure proposed by Zhao et aL^S 
but a small component from the i£f type Ce magnetic 
order induced by the Fe sublattice has to be included to 
fully explain our data below T^ e . 



2. PrFeAsO 

The ground state 3 H4 multiplct (J = 4) of the non- 
Kramers ion Pr 3+ in the C2 V low symmetry crystal field 
splits into nine singlets. Since the exact structure of the 
CEF levels is not known, the Pr 3+ magnetic susceptibil- 
ity cannot be estimated. The temperature dependence 
of the /xSR frequency in PrFeAsO is similar to the one 
observed in LaFeAsO. This suggests that the induced Pr 
moment is not as strong as in CeFeAsO. However, the 
comparison of the high temperature fit of the /iSR fre- 
quency for La, Ce, Sm, and Pr compounds yields the 
lowest value of /o = 22 MHz for the Pr compound. Ad- 
ditionally, there is a small but systematic deviation be- 
tween ^tSR and Mossbauer data for PrFeAsO in the Fe 
AFM phase. This demonstrates that also in PrFeAsO a 
Fe induced ordered Pr moment is present in the Fe AFM 
phase. One can conclude that, in contrast to the Ce 
compound, the Fe-Pr coupling constant I^~ Pr (Ku) is 
small and negative. According to equation the fields 
created by the Fe order parameter L^ x and the induced 
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Pr order parameter L£, r then have to point in opposite 
directions at the muon site. Therefore the /iSR frequency 
in PrFcAsO in the Fe AFM phase is smaller compared to 
the one in the other compounds. 

Recent neutron diffraction studies reveal that in the 
Pr ordered phase the Pr moments are directed along the 
c-axis, which creates a Lfj non-zero magnetic order pa- 
rameter, i.e. the same basis function of the irreducible 
representation as the induced Pr moments in the Fe or- 
dered phased Therefore, it is reasonable to assume that 
also in the Pr-ordered state, the local field at the muon 
site created by the Fe and the Pr magnetic order have op- 
posite signs. Respective magnetic structures are shown 
in Fig. [5Jb). The magnetic structure shown can there- 
fore qualitatively explain the 2 MHz drop of the /zSR 
frequency in the Pr ordered phase. However, the quanti- 
tative evaluation of the fiSR frequency using © and the 
Fe and Pr moment of 0.48 /ze/Fe and 0.84 /Us/Pr at 5 K 
that were deduced from neutron measurements^ reveal 
a strong discrepancy. If the suggested magnetic structure 
would be realized a /iSR frequency near to zero should 
have been observed. Even if one includes a strong in- 
crease of the local Fe moment (neglecting all transferred 
hyperfine field contributions) below the Pr ordering tem- 
perature which might be suggested by the Mossbaucr 
data^ the calculated /iSR frequency is much too small. 
To remove the discrepancy between neutron and /iSR ob- 
servations one has to conclude that in the PrFcAsO the 
Pr magnetic structure is either easy-plane non-collincar 
or collincar. In this case the situation would be simi- 
lar to the Ce and Nd compounds where in addition to 
the L^l{K.ii) order parameter also other a6-planc com- 
ponents like -Lf^T(Krj) or L^y(Kjj) exist. Similar con- 
clusions have been drawn from recent neutron diffraction 
studies*^ The consequences of non-collinear order pa- 
rameters will be discussed in Sec. I VIII 



3. SmFeAsO 

Compared to La and Pr, the Sm compound demon- 
strates a very similar /iSR response in the Fe ordered 
AFM phase which implies a L3 X (Kj) or L^CKjj) 
magnetic Fe order parameter. In spite of the similarity 
of the lowest CEF levels in Sm and Ce compounds the 
low value of the Sm g-factor reduces the induced Sm 
magnetic moment and therefore its contribution to the 
/iSR frequency in the Fe AFM phase. While in the Sm 
ordered phase, in contrast to the Pr and Ce compounds, 
the fiSK spectra change drastically and several well 
resolved frequencies are observed. At the lowest tem- 
peratures three /iSR frequencies are observed: 15 MHz, 
23 MHz, 31 MHz for which the approximate ratio of 
amplitudes is 1:4:1. This behavior indicates that for 
the given magnetic symmetry of the Sm order, at least 
three muon sites are inequivalcnt. This inequivalence 
can be caused by different translational symmetry of 
Sm and Fe magnetic order parameters or a complex 



non-collinear magnetic order of the Sm subsystem. 
A minimal model which can explain the /xSR spectra 
involves a L^ x (Ki) order parameter for the Fe sub- 
system with Ki and 1/2 [Lf™(K n ) +L§™(K n )] and 
1/V8 [Lf™(K n ) - Lf-(K„) + Lf-(K„) - Lf ?(K n )] 
order parameters for the Sm subsystem with Kjj trans- 
lational symmetries. This magnetic structure is shown 
in Fig. [UJc). It results in three inequivalent muon sites. 
The first with a low local field at the muon positions 
(15 MHz at A5 and A7), the second where the field from 
the R system cancels (23 MHz at A2, A4, A6, and A8) 
and a third with a high local field (31 MHz at Al an 
A3). The experimentally observed amplitudes are also 
reasonably well reproduced by the model. Within this 
magnetic model the size of the Sm magnetic moment 
can be estimated to be 0.4 /Lts/Sm at 1.9 K. Note that 
it is possible to fit the low temperature data with even 
more frequencies. The present statistical accuracy is not 
sufficient to determine the correctness of such a fit, but 
in a recent /iSR study the low temperature data have 
been fitted with five frequencies.— One has to note that 
our symmetry analysis also allows more complicated and 
even incommensurate Sm magnetic structures^ 

VII. RARE EARTH NON-COLLINEAR ORDER 
IN iiFEASO 

The above proposed R magnetic structures for the Ce, 
Sm and probably also in the Pr compound involve the 
well known exchange non-collinearity, i.e. a perpendic- 
ular orientation of neighboring R moments. As shown 
above, such a magnetic structure can be described by a 
composition of the Cartesian components of two different 
permutation modes like L^ x (Kuj)+Lf y (Kij) in the case 
of the Ce compound or like 1/2 [Lf z (Kj/) + Lg(K u )] 
in the case of SmFeAsO. From a pure symmetry point 
of view, commensurate and non-collinear magnetic or- 
der is characteristic for crystallographic structures with 
higher than orthorhombic symmetry that implies exis- 
tence of two- or three-dimensional IR for their space 
groups at k=0. From a pure energy point of view, the 
exchange non-collinearity results in competition of two 
or three different type of exchange multiplets (permuta- 
tion modes) which belong to the same IR and therefore 
possess the same Heisenberg energy. In high symmetry 
magnets this competition is resolved by accounting for 
fourth or higher order magnetic interactions in the sys- 
tem's free energy 45 i 46 Due to its crystallographic struc- 
ture the condition for accidental energy degeneracy of 
L R and L% rare earth permutation modes is fulfilled in 
the i?FeAsO compounds. The part of the Hamiltonian 
containing the R-R exchange interactions for Kn trans- 
lational symmetry can be written in the form: 

H?- R = ... + J«(L«f + J«(L§f + ... (15) 
where 
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J? = K R ~ R 
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For simplicity we omit the notation Kji. The symbol 
K R ^~ R denotes an exchange interaction coupling constant 
between nearest neighbor rare earth ions on sites 1 and 
(see Fig. [11] To deduce equation (fTB")) we used a permu- 
tation symmetry in which the interaction between ion 2 
and 4 is equal to the interaction between 1 and 3 and so 
on. If the exchange constant J R is negative and less than 
the others J R (j3 ^ a) a L R type of exchange magnetic 
order is realized, i.e. the relative orientation of magnetic 
moments (parallel/ anti parallel) in the magnetic struc- 
ture yields a non-zero Lf order parameter as given by 
equation ([6]). 

Neutron diffraction reveals a collinear if -type of rare 
earth exchange magnetic order in PrFeAsO non- 
zero magnetic order parameter)^ and in NdFeAsO (su- 
perposition of L^f and L^ z d non-zero magnetic order 
parameters) 4£ Similar magnetic structures have been 
proposed by a theoretical investigation of magnetic or- 
der in i?FeAsOfor R = Ce, Pr4£ The authors of Ref. \M 
claim that the R magnetic ground state is composed from 
two adjacent zigzag rare earth chains running along the 
a-axis which carry alternating up and a down spins. In 
our notation this is the L R type of exchange magnetic 
order. However, for the Ce compound this theoretical 
conclusion is not compatible with experimental data. 

A possible reason for the observed non-collinearity is 
provided by structural features of i?FeAsO. Note that 
without orthorhombic distortion, i.e. in the tetrago- 
nal PA/nmm phase, two exchange constants and 
K R ^~ R are equal to each other, and therefore their con- 
tributions to J R and J R cancel. Moreover, even in the 
orthorhombic Cmma phase their impact should be small. 
Indeed, the sign of R-R nearest-neighbor (NN) super- 
exchange depends on the R-O-R angle. There is a crit- 
ical value of this angle at which the exchange coupling 
constant changes its sign. For instance, for Cu-O-Cu 
bond this critical angle is equal to 104° 42 The exchange 
interaction is antiferromagnetic for magnetic bonds with 
angles higher than this critical angle. In the Ce com- 
pound the angles of Ce(l)-0-Ce(2) and Ce(l)-0-Ce(4) 
magnetic bonds, depicted in the Fig. [TT] as a 2 and a^, 
have values of a 2 = 105.39° and a A = 105.85° at T = 
1.4 K while the Ce(l)-0-Ce(3) angle is a 3 = 117.49°4£ 
For NdFeAsO the respective angles are: a 2 = 105.13°, 
a 4 = 105.58°, and a 3 = 118.06° at T = 0.3 K4I 

The angles a 2 and are always very similar in 
i?FeAsO. This implies that the exchange constants are 
also similar even if the value is far from the critical an- 
gle. The value is far from the critical angle that has been 
determined for the Cu-O-Cu bond, but might well be 




FIG. 11: The angle enumeration for R-O-R magnetic bonds 
in the rare earth tetrahedron of iJFeAsO (R = Pr, Sm, Ce). 



near critical for the R-O-R magnetic bonds considered 
here. The difference between Ce(l)-Ce(2) and Ce(l)- 
Ce(4) distances is only 0.011 Therefore, mainly anti- 
ferromagnetic exchange K R 3 ~~ R > determines the mag- 
nitude of the J R and J R coupling constants. In the limit 
of J R = J R — J R the corresponding exchange Hamilto- 



n R ~ R = j r [(Lf f 



+ (L?) 2 ] 



= 1/8 J R [(mf - mf ) 2 + (m R - mf ) 2 ](18) 

Using equation (|18| one can qualitatively explain the 
features of R magnetic order in the i?FeAsO compounds. 
In the R-O-R layer both planes of the R ions order anti- 
ferromagnetically with non-zero antiferromagnetic vec- 
tors if = mf — mf and if = mf — mf for R- 
planes which are located above and below the oxygen 
plane respectively. However, these two rare earth planes 
are not magnetically coupled, i.e. the mutual orienta- 
tion of those vectors is not defined by a Hcisenberg ex- 
change interaction. This magnetic frustration can be 
lifted by higher order interactions like biquadratic ex- 
change D R [(L R ) 2 - (Lf ) 2 ] 2 cx D R (l R ■ If ) 2 4£ For the 
case of J R < < D this results in a non-collinear mag- 
netic ground state with lf _Llf . Note that this is exactly 
the spin structure proposed for SmFcAsO in this study, 
and that it is very similar to the one proposed by Zhao et 
al4£ for CeFeAsO which is also consistent with our data. 
For the case J R < 0, D < a collinear structure with 
lf || lf is obtained, this structure is realized in NdFeAsO 
and probably also in PrFeAsO. The actual orientation of 
the R magnetic moments is determined by the complex 
hierarchy of second order anisotropy and fourth order 
interactions^ The orthorhombic distortions observed in 
the -RFeAsO system should lead to a weak coupling of 
antiferromagnetic R planes below and above the oxygen 
plane. Up to now only the R-O-R exchange paths have 
been considered. Note that an inclusion of the R-As-R 
exchange paths does not change the proposed scenario 
because the R-As-R magnetic bonds only contribute to 
the K^f 11 exchange and not to the R(l)-R(2) and R{1)- 



i?(4) exchanges. The hypothetical case where the ferro- 
magnetic i?(l)-As-i?(3) exchange cancels the antiferro- 
magnetic R(l)-0-R(3) exchange is very unlikely. 

VIII. INFLUENCE OF THE R MAGNETIC 
ORDER ON THE FE SUBLATTICE 

In the above discussion we focused on the action of the 
Fe order on the R subsystem. In accordance with Table 
HI! and the Hamiltonian ©, the L 3x Fe order parameter 
creates a molecular held on the R sites which is directed 
along z-axis and has a L R Z symmetry. The reverse ac- 
tion of the R magnetic order on the Fe subsystem should 
also be taken into account. This action will be especially 
important for a realization of the easy plane anisotropy 
in the R subsystem when the rare earth and iron order 
parameters possess different symmetry. The onset of the 
L R X or L R y in plane R order parameters will induce an 
out of plane canting of the iron moments along the z-axis 
due to these Fc-i? interactions: 
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FIG. 12: Different types of canted Fe structures that are in- 
duced by different rare earth magnetic structures, (a) Out of 
plane non-collinear Fe order {L$ x cos9 + Li z sin 9) induced by 
Z/f H -type rare earth order; (b) out of plane collinear Fe order 
(Lzx cos 9 + La? sin#) induced by L^-type rare earth order; 
(c) in plane non-collinear Fe order (Ls x cos 4> + L\ y sin <j>) in- 
duced by Lf z -type rare earth order; (d) in plane collinear Fe 
order (Lz x cos <j> + Ls y sin </>). Here, 9 and (j> are the polar an- 
gles between the Fe(l) magnetic moment and the z-axis and 
the x-axis, respectively. 



Kan-fx = ■ ■ ■ + l[xz i? ( K //) i 3x(K//)L| f 2 (K7/) + 

+I^- R (K II )L 3z (K II )L*(K II ) + 
+l[;- R (K II )L u (K II )Lf y (K II ) + 
+l[;- R (K II )L ly (K II )L R z (K II ) + ... 

(19) 

Interestingly, on the one hand the non-collinearity in 
the R system (L R Z and L R y order parameter) also induces 
a non-collinearity in the Fe order (L 3x and L\ z order pa- 
rameters). On the other hand, a collinear and oblique R 
magnetic order (L R X and L R Z order parameters) induces 
a Fe magnetic structure which is still collinear but the Fc 
moments deviate from the x-direction. Different types 
of canted Fe structures are shown in Fig. [T^l 

We discuss now whether a canting of the magnetic 
structures could be observed by /iSR and Mossbauer 
spectroscopy. Here we consider the contribution from the 
Fe sublattice only. Supposing that the size of the Fe mag- 
netic moment does not change below the R ordering tem- 
perature, i.e. it is close to saturation, the /iSR response 
is affected by the canting in the following way: From the 
dipole held calculation presented in Appendix A, it fol- 
lows that the magnetic held at the muon site created by 
both non-collinear Fig. [T27a) and collinear Fig. [T27b) 
out-of-plane structures is tilted by the out-of-plane an- 
gle 6 but does not change its modulus. Therefore, ^SR 
experiments on powder samples cannot distinguish the 
canted structures from the original collinear L 3x struc- 
ture. In contrast, the held at the muon site, i.e. the ^SK 
frequencies depend on the angle <f> for in-plane canted 
structures shown in Fig. [T27 c) and (d). In the case of 
the magnetic structure shown in Fig. [TSTd) the /iSR fre- 
quency always decreases like cos <fi, whereas, in the case of 
structure (c) in Fig. [12] the field at the muon site may de- 
crease or increase depending on the sign of <f> and has the 



following angle dependence V2 cos((/> + 7r/4). In turn, its 
sign determines the sign of the coupling constant I^y^ 11 
of the Fe-R interaction. 

Mossbauer spectroscopy, in contrast to yuSR, should 
be very sensitive to the orientation of the Fe moments. 
Our calculation of the electric field gradient tensor on 
Fe site reveals a strong anisotropy with the largest com- 
ponent in the z-direction. This makes Mossbauer spec- 
troscopy a suitable tool to detect the direction of the 
Fe magnetic moment especially in the case where out 
of a6-plane deviations from x-axis are present. Again 
supposing an almost constant modulus of the iron mo- 
ments below T R one can derive a strong increase of the 
quadrupole shift along with an increase of the average 
hyperfine field for the out of plane canted iron structures 
like those shown in Fig. [T2Ta') and (b). This behav- 
ior has already been detected in the Mossbauer studies 
of PrFeAsO^ McGuire et al. come to the conclusion 
that the Fe spin-reorientation with a significant compo- 
nent of Fe moment along the c-axis takes place below the 
Pr ordering temperature. For the in-plane collinear and 
non-collinear canted structures shown in Fig. [TSTc) and 
(d) the dependence of Mossbauer spectra on the rotation 
angle (f> should be weaker than for out-of-plane rotation 
since it only changes the asymmetry parameter of the 
electric field gradient tensor. 

Combining both comments above, one can understand 
why the drastic change of PrFeAsO Mossbauer spectra 
below T^ r has no analogy in the /iSR spectra. The fol- 
lowing model can unify both observations: 

i) In the Pr subsystem, the L R type of exchange order 
is realized with an easy plane anisotropy. 

ii) The order parameter which appears below T^ r is 
L R X with a small induced L R Z component remaining even 
down to the lowest temperatures due to the Fe molecular 
field I R x e ~ Pr (Kii)L 3x (Kn). The latter component being 
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negative as mentioned above. 

iii) Below T^ r the in-plane collinear Fe magnetic struc- 
ture gradually rotates into the collinear out-of-planc 
L^x cos 9 + L^ z sin 9 structure. 

Respective magnetic fields at the muon sites expressed in 
MHz are given by 

B X {A) = f sin 9- 13.55- 0.59- L^(Ku) 
B Z {A) = / o cos0. (20) 

Here /o is the low temperature frequency extrapolated 
from temperatures well above T^ r in the Fe antifcrro- 
magnetic phase. It is determined by the saturated Fe 
magnetic moment and has a value of /o = 22 MHz in 
the Pr compound. In (|20p the small contribution from 
the L^ z component which decreases as the ordered iron 
moment rotates is omitted. The decrease of the /xSR 
frequency below T^ r is roughly 2 MHz. Assuming a 
magnetic moment of approximately Iiib/Fy a canting 
of sin# = 0.42 is obtained from Eq. (|2"0|), This esti- 
mation confirms indeed a strong tilting of the Fe mo- 
ments to the z-axis which was already concluded from 
Mossbauer spectroscopy experiments^ The respective 
magnetic structure is shown in Fig. (JSJd). Note, that 
in this model the angle 9 can achieve 90°, if we sup- 
pose the Pr moments equal to 0.26/is/Pr. As mentioned 
above, the iiSFL signal would be exactly the same for an- 
other combined Fe-Pr magnetic structure with a non- 
zero non-collinear L^ x cos 9 + L\ z sin9 in the Fe and an 
Lfy order parameters in the Pr subsystem. We conclude 
that these two different magnetic structures are indistin- 
guishable by fiSR and Mossbauer studies. 

Finally, we comment on other aspects of the Fe-R in- 
teractions. It is generally known that the Fe-R inter- 
action plays a decisive role for the magnetic anisotropy 
of the Fe subsystem in Fe-i? compounds. For instance, 
temperature induced spin reorientation phase transitions, 
even far above of rare earth magnetic order temperature, 
have been observed in rare earth orthofcrritcs ReFoOa 51 
and in the ReFenTi compounds.— The temperature de- 
pendence of the Fe magnetic anisotropy is determined by 
competing contributions from the Fe and R subsystems. 
In ReFeC>3 compounds the R ordering is almost always 
accompanied by a Fe spin reorientation, and the canting 
angles are relatively large due to a strong rcnormaliza- 
tion of the Fe magnetic anisotropy. In this sense it would 
be surprising if the i?FcAsO system were an exception 
to this rule. Furthermore it is clear that the large pro- 
posed Fe canting angles 9 in Pr ordered phase of PrFeAsO 
can be achieved by a renormalization of the Fe magnetic 
anisotropy but not solely by the direct Fe-Pr interaction 
mechanism (fl9")) . 

Note that our fj.SK studies of i?FeAsO with Sm, Pr 
and Ce do not show temperature induced spin reorienta- 
tions of the iron magnetic structure in the Fe AFM phase. 
Spin reorientation starts at a temperature at which the 
anisotropy changes sign, i.e. it develops as a standard 



second order phase transition. These have not been de- 
tected, neither by change of the muon spin relaxation 
rate nor by other methods like specific heat or suscep- 
tibility studies. In the low temperature range the iron 
spin reorientation begins with the onset of the rare earth 
magnetic order and therefore it can be masked by the 
latter. 



IX. CONCLUSION 

In conclusion we presented a detailed study of the mag- 
netic order and the interplay of the rare earth and iron 
magnetism in i?FcAsO with R = La, Ce, Pr, and Sm, 
the magnetically ordered parent compounds of the re- 
cently discovered i^Oi-^F^FeAs high temperature super- 
conductors. Using zero field (ZF) /iSR, the Neel temper- 
atures, as well as the temperature dependence of the sub- 
lattice magnetizations, have been determined with high 
precision. 

The second order Fe magnetic phase transition is well- 
separated from a structural transition which occurs 10- 
20 K above TV. The muon site in the i?FeAsO crys- 
tal structure has been obtained by electronic poten- 
tial calculations involving a modified Thomas-Fermi ap- 
proach. Using calculated dipolc fields at the muon site, 
our /iSR experiments indicate an antiferromagnetic com- 
mensurate order with the iron magnetic moments di- 
rected along the a- axis above T§. This is consistent 
with neutron scattering results reported earlier^ 3 - The 
calculations show, that the two types of iron order pa- 
rameter with different translational symmetry, that have 
been observed in i?FcAsO compounds, i.e. with (K/) 
and without (K-n) doubling of the magnetic primitive 
cell along the c-direction cause nearly the same field at 
the muon site and are therefore not distinguishable by 
^SR. Mossbauer spectroscop y 1 ^ 23 as well as fiSR prove 
that the ordered Fe magnetic moment is approximately 
0.4 fiB and does not vary by more than 15 % within the 
series i?FcAsO with R = La, Ce, Pr, and Sm. This is in 
sharp contrast to published neutron results^ 3 - where Fe 
magnetic moments of 0.36 fiB, 0.48 /ib, and 0.8 fiB have 
been deduced for R = La, Pr, and Ce respectively. In 
the neutron studies, the observed intensity in the mag- 
netic Bragg peaks is solely attributed to an ordered Fe 
moment for temperatures larger than Tjy . Here we could 
show that a sizable magnetic moment of the R subsys- 
tem which is induced by the ordered Fe subsystem is 
already detectable far above T$. Especially in CeFeAsO 
this R magnetization is exceptionally large and amounts 
to approximately 0.3 ns/Ce just above T^ e . Since for 
symmetry reasons the Ce and Fe magnetic moments con- 
tribute to the same Bragg peaks, neglecting the Ce mag- 
netization results in the strong overestimation of the Fe 
moment in the CeFeAsO system by neutron scattering. 

Static magnetic ordering of the rare earth moments has 
been observed for the R = Ce, Pr, and Sm compounds 
with transition temperatures of Tjy = 4.4(3), 11(1), and 
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4.66(5) K respectively. Using available literature data of 
Mossbauer spectroscopy and neutron diffraction and the 
[iSR data presented here, we propose combined R and Fe 
magnetic structures below Tjy for all investigated com- 
pounds. For CeFeAsO, the //SR data are consistent with 
a non-collinear easy plane AFM order of the Ce moments 
as concluded from neutron diffraction.— For SmFeAsO, 
several new /iSR frequencies develop in the Sm ordered 
phase. To explain the /^SR spectra a minimal model of 
non-collinear Sm magnetic order is proposed. The ad- 
ditional /iSR frequencies originate from different trans- 
lational symmetries of the Sm and Fe order parameters. 
Analyzing structural features of i?FeAsO we argue that 
the non-collinear rare earth magnetic order observed in 
Ce and Sm compounds arises due to a weak magnetic 
coupling of the adjacent R planes in the R-O-R layer. 
A weak coupling of the R-O-R layers is probably an in- 
herent feature of all i?FeAsO compounds. For PrFcAsO, 
we propose that the Pr moments order also with an easy 
plane anisotropy. Our model therefore differs from the 
one determined by Zhao et al.— which implies an easy 
axis anisotropy in the Pr ordered phase. Only the model 
with easy plane anisotropy is able to consistently describe 
the /xSR and Mossbauer— as well as the neutron data. 

Even though Fe spin reorientation phase transitions are 
frequently observed in R— Fe systems, e.g. in rare earth 
orthoferrites RFeO^ or i?FenTi^, we do not find such 
a transition in i?FeAsO above T$ for R = Ce, Pr and Sm. 
In the R ordered phase one can expect a strong influence 
of the rare earth magnetic subsystem on the Fe ordering 
direction especially in the case when the iron and rare 
earth order parameters have different symmetry. This is 
the case for collinear as well as non-collinear easy plane 
AFM order in the R subsystem. We suggest that this is 
the reason for the Fe spin reorientation developing below 
T^ r as deduced from the /LtSR data on PrFeAsO. 

Using symmetry arguments we demonstrate the ab- 
sence of a Heisenberg magnetic interaction between the 
Fe and the R subsystem in i?FcAsO. Therefore, the ap- 
parent Fc-i? interaction is realized by a non-Hciscnbcrg 
anisotropic exchange. Our calculations additionally show 
that dipole-dipole interactions are much too weak to ac- 
count for the observed couplings. We showed that in Ce- 
FeAsO the Fe-Ce coupling is exceptionally large. This to- 
gether with the large paramagnetic Ce magnetic moment 
explains the sizable Ce magnetization observed in our 
experiment in the Fe AFM ordered phase as well as the 
presence of a small component of Ce moment along the c- 
axis in the Ce AFM ordered phase. For PrFeAsO, from 
a comparison of /xSR, Mossbauer and neutron diffrac- 
tion data we derive a much weaker but noticeable Fe-Pr 
coupling constant which has opposite sign compared to 



the Fe-Ce coupling constant. For SmFeAsO the Fe-Sm 
coupling is very weak and the induced polarization of 
the Sm moments by the magnetically ordered Fe sys- 
tem is not measurable above T^ m . Even though the 
lowest CEF levels in Sm and Ce are relatively similar, 
the low value of the Sm g-factor reduces the induced Sm 
magnetic moment and therefore its contribution to the 
/iSR frequency. The exceptionally strong coupling of the 
Ce to the Fe subsystem in CeFeAsO can also be under- 
stood from band structure calculations: Only CeFeAsO 
possesses a considerable 3d-4f hybridization . 24 ! 25 ' 26 Ad- 
ditionally, the strong interaction of the two magnetic sub- 
lattices is in accordance with the observed strong electron 
correlation of the Ce 4f electrons and the heavy fermion 
behavior observed in the related phosphide CeOFeP.— 

Finally, from our analysis we can conclude that the 
magnetic i?-Fe interaction is probably not crucial for the 
observed enhanced superconductivity in ROi-xF x ¥eAs 
with magnetic 4f ions compared to LaOi-^F^FeAs, since 
only in CeFeAsO a strong R-Fe coupling is observed, 
while it is much weaker in the Pr and Sm compounds. 
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X. APPENDIX 

A. Magnitude and symmetry of dipole fields 
created by the iron and rare earth subsystems at the 
A type muon site 

Below we provide the calculation of magnetic dipole 
fields at the muon site A for different magnetic modes and 
different translational symmetries. The magnetic fields 
below are given in units of 10 _1 T and basis functions L 
in units of hb- The coordinates of the muon site A are 
(0,1/4,0.41) 



19 



B X (A) 

By (A) 

B Z {A) 

-0.79 
-0.66 





-2.56 
-2.56 





2.58 \ /L te (Kj) 
L 3y (K!) 
2.58 / V L 3z (K T ) 



Lg(K T ) 




B X {A) \ 

By (A) 

BAA) J 



-0.89 






/ L lx (K H 
V L lz (K H 




-0.59 
1.41 



fL? x (K H ) 




2.57 \ / L 3x (K u ) 
L 3v (K u ) 
2.57 / I L 3z (K H ) 




fLg(K H ) 
L§y(Kjj) 
\Lg(K H ) 



B X (A) 




By (A) 


)-(? 


B Z (A) 




-0.40 








-0.40 





0.80 



2.12 



f L lx (K 



£fx(K///) 




III) \ 







in) 


*(■ 


-4.44 


hi) J 







-0.43 











-0.43 











0.86 



-4.44 \ / L 3x (K ln ) 
L 3y {K IH ) 
J \L 3z {K IH ) 



B X {A) 
By (A) 
B Z {A) 

-0.95 







-0.95 
1.90 



^(K ) 
F Z R (K ) 



F X (K Q ) 
^(Ko) 
i^(Ko) 



-4.45 \ / L 2x (K ) 
-4.45 L 2y (K ) 
/ \ L 2z {K ) 



-0.44 
-0.44 
0.88 



i|(Ko) 




r 



B. Dipole fields created by the iron and rare earth 
subsystems at the rare earth site 

Below we provide the calculation of magnetic dipole 
fields at the rare earth site for different magnetic modes 



and different translational symmetries. The magnetic 
fields below are given in units of 10 _1 T and basis func- 
tions L in units of ^b- 
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